Modeling of Fe pnictides: the Magnetic Order and Pairing Channels 
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A two-orbital model for Fe-pnictide superconductors is investigated using computational tech- 
niques on two-dimensional square clusters. The hopping amplitudes are derived from orbital overlap 
integrals, or by band structure fits, and the spin frustrating effect of the plaquette-diagonal Fe-Fe 
hopping is remarked. A spin "striped" state is stable in a broad range of couplings in the undoped 
regime, in agreement with neutron scattering. Adding two electrons to the undoped ground state of 
a small cluster, the dominant pairing operators are found. Depending on parameters, two pairing 
operators were identified: they involve intei-xz-yz orbital combinations forming spin singlets or 
triplets, transforming according to the _B2g and A2g representations of the _D4h group, respectively. 



Introduction: The recent discovery of superconduc- 
tivity in the layered rare-earth oxypnictides compounds 
LnOj_3;Fa;FeAs (Ln=La, Pr, Ce, Sm) has captured 
the attention of the condensed matter community 
The high current record critical temperature Tc ~55 K 
in SmOi_2;F3;FeAs [§] suggests that an unconventional 
pairing mechanism may be at work 3, 4]. 

As for Cu-based high temperature superconductors 
(HTSC), the analysis of undoped compounds, such as 
LaOFeAs, is expected to provide important informa- 
tion toward the understanding of the superconducting 
(SC) state reached by ^10% F doping. Neutron scat- 
tering experiments have provided evidence of magnetic 
order in LaOFeAs at 134 K: Fe spins order into fer- 
romagnetic "stripes" that are aligned antiferromagneti- 
cally 5-7]. In the two-dimensional (2D) square lattice 
notation, the LaOFeAs magnetic structure factor has 
peaks at wavevectors 9^(0, 7r),(7r, 0) Assuming a 

smooth continuity between the undoped and F-doped 
compounds, the pairing mechanism could be magnetic 
in origin and triggered by this unusual magnetic state. 

Theoretical work on the new superconductors includes 
band structure calculations that have shown the rele- 
vance of the 3c? levels of Fe [3, . A metallic state involv- 
ing a Fermi surface (FS) made out of disconnected small 
pieces ( "pockets" ) was predicted [sj . To understand some 
of the properties of the undoped limit, electron correla- 
tions appear to be important [loj . Two-orbitals descrip- 



tions [Hi, and other models have been proposed, and 
a variety of approximations have lead to several uncon- 
ventional pairing channel proposals 13- 15|. 

Our purpose is to report the first unbiased computa- 
tional results obtained using a model Hamiltonian for Fe 
pnictides, with emphasis on a real-space description. In 
the undoped limit, a spin striped magnetic state is ob- 
tained and explained. With light electron doping, novel 
pairing operators are identified. The path followed here 
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FIG. 1; (Color online) (a) Small cluster illustrating the ge- 
ometry of the FeAs layer. Open (filled) green circles indicate 
As positions above (below) the Fe plane, (b) The Fe-Fe NN 
path, (c) The Fe-Fe NNN path, (d) Fermi surface of //foAs for 
pdn/pda=— 0.2 in the U = J = limit, (e) Same as (d) but 
for the parameters of Ref. [111]. In (d), the half- filled chemical 
potential is at —0.03. (f) The tilted 8-site cluster used here. 



mimics research in the HTSC, where the computational 
study of model Hamiltonians in real space [l6j provided a 
perspective dual to momentum-space diagrammatic cal- 
culations. In fact, early numerical studies of the 2-hole 
state on small t-J clusters indicated that the pairing was 
in the d^2_y2 channel [l6| . Thus, it is natural to follow a 
similar path for the new Fe superconductors. 

Model and Techniques: Fe pnictides have a layered 
structure, with the Fe atoms forming a 2D square lattice 
and the As atoms located above or below the plane, at the 
centers of alternating plaquettes, see Fig. 1(a). Here, the 
emphasis will be on the Fe dxz and dyz degenerate states 



111 Il2j since band structure calculations have shown the 
importance of these orbitals at the Fermi surface 0, 113] ■ 
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Since the complexity of the problem rapidly increases 
with the number of orbitals, it is reasonable to start with 
just two orbitals, contrast the results with experiments, 
and slowly build up a more realistic model with extra or- 
bitals. To estimate the hopping amplitudes, the Slater- 
Koster (SK) tight-binding scheme is followed [l^. This 
approach is simple, analytical, and leads to a geometrical 
understanding of the magnetic phase. The SK method 
for the hopping integrals needs as input the location of 
the Fe and As atoms, and the nature of the orbitals. The 
Fe-Fe (Fe-As) distance used is r=2.854 A (s=2.327 A). 
The effective Fe-Fe hopping amplitudes - via As - are the 
product of the direct Fe-As hoppings, the three p orbitals 
were taken into account on the intermediate As ion. Two 
Fe-As-Fe paths connect nearest-neighbor (NN) Fe-sites, 
while only one exists for next-nearest-neighbor (NNN) 
Fe's along the plaquette diagonals. The kinetic energy of 
the resulting model restricted to Fe sites is 



The on-site Coulombic terms include a Hubbard re- 
pulsion U for electrons with the same a, a repulsion U' 
for different a, a ferromagnetic Hund coupling J, and a 
pair- hopping term with strength J'~J [20i |: 

Hint ni,a,^ni,a,i + {W - J/2) ^ rii^a^rii^y 

i,a i 
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H. c. 



where ^ ^ creates an electron with spin cr in the 
orbitals a=x,y (dxz and dyz, respectively) at site i 
of a 2D square lattice. x and y are unit vectors 
along the axes. The SK-evaluated Fe-Fe hopping am- 
plitudes are ti=-2[{b'^ - a^) -t- g^l/Apd, t2=-2[{b'^ - 
a^) - g^]/Apd, h=-{a^ + ~ 9'^)/\d. and U^-{ab - 
g^)/Apd, where the Fe-As hopping amplitudes are 
a=0.324(pdcr) - 0.374(pd7r), 6=0.324(pdcr) + 0.123(pd7r), 
and g—0.26i{pda) + 0.31(pc?7r). pda and pdn are SK 
parameters and Apd is the energy difference between 
the p and d levels. The overall energy scale is set by 
{pda)'^ / Apd, which is of the order of eV and will be used 
as unit of energy, pdn/pda is a free parameter in . 
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Equation (1) is formally the same as in Refs. 
but the values for the hoppings are different: our ap- 
proach relies on the analytic calculation of the hoppings 
in a "first-principles" SK-based context, while Ref. |11| 
fits the hoppings to bands from LDA calculations. An 
interesting conclusion of our effort is that both sets of 
parameters lead to similar results for the magnetic order 
and the pairing. Eq uation (1) has invariance under the 
D4h point-group [l9|, including a tt/2 rotation of the lat- 
tice together with orbital exchanges a; — > y and y — > —x. 



Si.Q (^^i.a) is the spin (density) in orbital a at site 
i. The standard relation U'=U — 2J due to rota- 
tional invariance was used [2lj. The full model becomes 
HFcAs—Hii+Hint- Since the NN Fe-As-Fe bond angle 
^NN [Fig. 1(b)] is closer to 90° than the NNN Fe-As- 
Fe angle 6'nnn [Fig. 1(c)], we find a strong NNN hopping 
ts and the ratio t^/ti is of order 1 for broad ranges of 
pdn/pda, without fine tuning. At intermediate to large 
U, the resulting effective Fe-Fe spin interaction along the 
plaquette diagonals consequently becomes as large as be- 
tween NN Fe sites [13], or even larger. In the early days 
of HTSC, investigations of the resulting frustrated effec- 
tive spin model unveiled a spin striped phase in the one- 
orbital model [13]. As shown in Figs. l(d,e), the non- 
interacting system has electron FS around the {X,Y) 
points, and hole FS around the F and M points, which 
are equivalent on folding the Brillouin zone. 2-orbital 
models cannot have both hole pockets around F point, 
as found in band structure calculations [1, [lH . 

Magnetic properties in the undoped limit: We study 
the ground state of model Eqs.(l,2) in the undoped hmit 
by using two techniques: Exact Diagonalization (ED) 
and the Variational Cluster Approach (VGA). The first 
method allows for an unbiased analysis, albeit restricted 
to small clusters [13] , while the second extends the calcu- 
lation to the bulk self-consistently [13, [2^ . We apply ED 
to the 2x2 and tilted VS x VS [Fig. 1(f)] clusters with 
periodic boundary conditions [3| ■ While the 2x2 cluster 
only has 4,900 states even if no symmetries are used, the 
8-sites cluster has 20,706,468 states with translational in- 
variance implemented, and is computationally demand- 
ing. We therefore fixed the ratio U/J to 4, compatible 
with some estimates However, other U/J rations (to 
be discussed in future publications) do not critically af- 
fect the results presented below. In particular, the spin 
striped state and the singlet pairing (discussed later) sur- 
vive for small J. The typical inequality \pdTT/pda\ < 1 is 
assumed, and the sign of pdn is chosen such that the FS 
agrees with band structure calculations (see below) . Both 
on the 2x2 and on the \/8 x ^/8, we observe magnetic or- 
der with q = (0,7r), (7r,0) in real-space spin correlations 
as well as in the magnetic structure factor S(q), see, e.g.. 
Fig. 2(a). This leads us to believe that size effects are not 
severe. The q^{0, tt), (tt, 0) state is stable at least in the 
large square —0.5 < pdn/pda < and < U < 4, and it 
is generated by the robust plaquette-diagonal hoppings. 
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FIG. 2: (Color online) (a) S{q) vs. g for the x clus- 
ter at the f/'s indicated, with J = 17/4 and pdir /pda=—Q.2. 
The momenta allowed in this cluster are (0, 0), (±7r/2, ±7r/2), 
(0, tt), (7r,0), and (vr, tt) [2H. (b) VGA grand potential [in en- 
ergy units (pd(T)^/Apd] vs. staggered magnetic fields ftstagg for 
q = (tt, tt) and (0, tt), and pd-K /pda— —0.2, U — 1, J = 0.25, 
The minimum for (0, tt) at /istagg 7^ indicates symmetry 
breaking, (c) mstagg vs. U, with U/J = 4. The small U 
region at pdn/pda=— 0.2 was numerically unstable. 
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FIG. 3: (Color online) One-particle VGA spectral function 
A(k,uj) of i^pcAs for (a) [7=0.5, J=0.125, and pdn/pda^-0.2, 
and (b) U—1, J=0.25, pd-K /pda— —0.5 in the symmetry bro- 
ken phase with (0, tt) magnetic ordering at half-filling. A 
broadening 0.05 was used, (c) Density-of-states obtained by 
k-integrating the spectral functions, at the U's indicated, with 
[7/ J=4 and pdvr /pd(T=— 0.2. (d) FS corresponding to case (a). 
The FS has been symmetrized under rotations. 



We find similar results for both the SK hoppings and 
those of Ref. [ll| , and our results also agree with weak- 
coupling RPA approximations, where similar order arises 
from nesting [ll|, [2^ . Further evidence that we indeed 
identified the dominant magnetic channel comes from the 
VGA [23] , where the self-energy of a small cluster is op- 
timized by varying appropriate "fictitious" fields such 
as chemical potentials or symmetry-breaking staggered 
magnetic fields 2J, |25| . It thus combines the exact solu- 
tion of a small cluster with access to the bulk limit [23| . 
The grand potential [Fig. 2(b)] demonstrates that the 
symmetry breaking indeed occurs in the q—(0,T:), (7r,0) 
channel. Figure 2(c) shows the stripe order-parameter 
"istagg (the staggered moment in units of Bohr magne- 
ton per Fe ion) vs. U: we find a (U, pdn/pda) regime 
that can accommodate the small mstagg found with neu- 
trons [5] . Note also that other experiments have reported 
larger Wstagg values [2]- 

Figure 3 shows photoemission spectra A{]i,uj). The 
first case (a) for pdn/pda = — 0.2, [/ = 0.5 has a dis- 
persion similar to that of the non-interacting system 
[Fig. 1(d)], its FS is shown in Fig. 3(d) and the nodal 
structure will be discussed in a future publication. The 
density of states (DOS) has a small pseudogap [Fig. 3(c)], 
somewhat deeper than for U = 0, while larger U = 2 
leads to an insulating hard gap. Another interesting 
regimes is shown in (b): At [/ = 1 and pdn /pdcr=— 0.5, 
the chemical potential lies in a region with many states. 



suggesting a correlated metal. 

Pairing channels: In the HTSC cuprates, the domi- 
nant pairing channel could be identified by adding two 
carriers to the undoped finite cluster (which has the same 
symmetries as the bulk system) and by then evaluating 
the quantum numbers under tt/2 rotations [l6j . Since 
NNN hoppings play a key role in the present case, we 
need at least a 2x2 cluster for each sublattice, the min- 
imal cluster satisfying the requirement is the \/8 x \/8 
cluster. Varying U and pdn/pda {U/J=4), we find that 
both singlet and triplet regimes can be reached by adding 
two electrons to the undoped (i.e. half- filled) system j28|], 
see the phase diagram in Fig. 4(a). At small \pdTT/pda\ 
and intermediate or large U, the total spin is 0, and in 
searching for the local operator connecting the ground 
states of the undoped and doped systems, the largest 
overlap at intermediate U is for 



(3) 



7t 



-a, -I' 



or in k-space, A'''(k)=^^(cos fca: + cosky)d]^ ^ .^d 
a = x,y and /t = x,y. This operator is a spin singlet that 
transforms as the representation of the D^oi group 
(lij . and it involves different x and y orbitals on NN 
sites to optimize the NN kinetic energy [Fig. 4(b)]. In 
other parts of the phase diagram, a spin-triplet domi- 
nates, which is odd under orbital exchange, transforms 
according to [l^ , and also involves different orbitals 
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FIG. 4: (Color online) (a) y/Ex y/8 cluster results with 2 more 
electrons than half-filling, at U/ J=4. Shown are regions with 
singlet and triplet pairing (see text). (b,c) Schematic repre- 
sentation of the singlet and triplet pairs (see text). Black ar- 
rows represent the magnetically ordered background. White 
arrows are the added electrons. x,y are the orbitals. 



on NN sites [Fig. 4(c)]. Its projection-1 operator is 



Ani)i = E(4.,t4 



(4) 



that in momentum space becomes A^(k)i = (cos fc^; + 

cos fcj,)(rfk,^,^c^Lk.,^,^-rfLy,t^-k,2;,t)- rescmblcs the op- 
erator of Ref. ,15] . although they use on-site pairing. Of 
the 16 possible pairing operators allowed by the symme- 
try of the Hamiltonian jiol . [29| . our singlet and triplet 
operators correspond to #9 and #12 of Ref. [19, 30]. 

Conclusions: We studied a simple model for the 
FeAs superconductors numerically. The undoped system 
shows q^(0, 7r),(7r, 0) spin order. We identified dominant 
pairing operators for two added electrons: depending on 
parameters they can be spin singlet or triplet, transform- 
ing non-trivially under ti/2 rotations. Future work will 
address more realistic models beyond two orbitals and 
contrast their results against those reported here. 

Acknowledgment: We acknowledge discussions with F. 
Reboredo. Research supported by the NSF grant DMR- 
0706020, the Div. of Materials Sciences and Eng., U.S. 
DOE under contract with UT-Batelle, LLC, and the Aus- 
trian Science Fund grant P18551-N16. 



[1] 

[2] 

[3] 

[4] 
[5] 
[6] 
[7] 

[8] 

[9] 
[10] 

[11] 
[12] 

[13] 



[14] 



[15] 

[16] 

[17] 
[18] 
[19] 
[20] 
[21] 
[22] 

[23] 
[24] 
[25] 
[26] 
[27] 

[28] 



[29] 
[30] 



Y. Kamihara et al, J. Am. Chem. Soc. 130, 3296 (2008). 

Z-A. Ren et al, Chin. Phys. Lett. 25, 2215 (2008); Eu- 

rophys. Lett., 83, 17002 (2008). 

A. Sefat et al, Phys. Rev. B 77, 174503 (2008). 

L. Boeri et al., Phys. Rev. Lett. 101, 026403 (2008). 

C. de la Cruz et al, Nature 453, 899 (2008). 

J. Dong et al, Europhys. Lett. 83, 27006 (2008). 

A.I. Goldman et al, Phys. Rev. B 78, 100506 (2008); C. 

Krellner et al, Phys. Rev. B 78, 100504 (2008). 

D. J. Singh and M.-H. Du, Phys. Rev. Lett. 100, 237003 
(2008); G. Xu et al, Europhys. Lett. 82, 67002 (2008). 
C. Cao et al, Phys. Rev. B 77, 220506 (2008). 

K. Haule et al, Phys. Rev. Lett. 100, 226402 (2008). 

S. Raghu et al., Phys. Rev. B 77, 220503 (2008). 

Q. Han et al., Europhys. Lett. 82, 37007 (2008); T. Li, 

J. Phys.: Condens. Matter 20, 425203 (2008). 

K. Kuroki et al., Phys. Rev. Lett. 101, 087004 (2008); 

I. Mazin et al., Phys. Rev. Lett. 101, 057003 (20 08); 

M.M. Kors hunov and I. Er emin, larXiv:0804.T793l G. 

Baskaran, arXiv:0804.1341' P. Lee and X-G. Wen, 

arXiv:0804.1739; Z-J. Yao et al, arXiv:0804.4166 , X-L. 

Qi et al., arXiv:0804.4332 

Q. Si and E. Abraham, Phys. Rev. Lett. 101, 076401 
(2008); C. Fang et al, Phys. Rev. B 77, 224509 (2008); 
and C. Xu et al, Phys. Rev. B 78, 020501 (2008). 
X. Dai et al., Phys. Rev. Lett. 101, 057008 (2008). See 
also W. Chen et al., arXiv:0808.3234 

E. Dagotto, Rev. Mod. Phys. 66, 763 (1994), and refer- 
ences therein. 

V. Vildosola et al, Phys. Rev. B 78, 064518 (2008). 
J. C. Slater and G. F. Koster, Phys. Rev. 94, 1498 (1954). 
Y. Wan and Q-H. Wang, arXiv:0805.09W, 
See e.g. A. M. Oles et al, Phys. Rev. B 72, 214431 (2005). 
Page 36 of E. Dagotto et al, Phys. Rep. 344, 1 (2001). 
C. Henley, Phys. Rev. Lett. 62, 2056 (1989); E. Dagotto 
and A. Moreo, Phys. Rev. Lett. 63, 2148 (1989). 
E. Dagotto et al., Phys. Rev. B 41, 9049 (1990). 
C. Dahnken et al., Phys. Rev. B 70, 245110 (2004). 
M. Potthoff, Eur. Phys. J. B 36, 335 (2003). 
T. Yildirim, Phys. Rev. Lett. 101, 057010 (2008). 
VGA is run on large but finite 48x48 lattices. We com- 
pared the (0, tt) phase to (tt, tt) and ferromagnetic order. 
As in "AF-broken-links" argumentations for HTSC [i6| , 
here pairing of carriers may emerge from the "damage" 
an added electron does to the magnetic background. 
Z-H. Wang et al, arXiv:08 05.0736i Y. Zhou et al., Phys. 
Rev. B. 78, 064514 (2008). 

In a future publication, we will show that similar results 
are also obtained using the parameters of Ref. [ij. 



